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ABSTRACT 

We report ground-based follow-up observations of the exceptional source, 
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ID141, one the brightest sources detected so far in the if- ATLAS cosmo logi- 
cal survey. ID 141 was observed using the IRAM 30-meter telescope and Plateau 
de Bure interferometer (PdBI), the Submillimeter Array (SMA) and the Ata- 
cama Pathfinder Experiment (APEX) submillimeter telescope to measure the 
dust continuum and emission lines of the main isotope of carbon monoxide 
and carbon ([Cl] and [ClI]). The detection of strong CO emission lines with 
the PdBI confirms that ID 141 is at high redshift {z = 4.243 ± 0.001). The 
strength of the continuum and emission lines suggests that ID 141 is gravita- 
tionally lensed. The width (AI^whm ~ 800kms~^) and asymmetric profiles 
of the CO and carbon lines indicate orbital motion in a disc or a merger. 
The properties derived for ID 141 are compatible with a ultraluminous (Lfir ~ 
8.5±0.3 X 10^^ /iL^ Lq, where fiL is the amplification factor), dense (n ^ 10^ cm~^) 
and warm (Tkin ~ 40 K) starburst galaxy, with an estimated star-formation 
rate of (0.7 to 1.7)xlO^/iL^ MQ/yr. The carbon emission lines indicate a dense 
(n Ri lO^cm"^) Photo-Dominated Region, illuminated by a far-UV radiation field 
a few thousand times more intense than that in our Galaxy. In conclusion, the 
physical properties of the high-^ galaxy ID 141 are remarkably similar to those 
of local ultraluminous infrared galaxies. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: high- 
redshift — galaxies: starburst — submillimeter: galaxies — individual (ID 141) 



1. Introduction 

Over the last decade, submillimeter surveys have revolutionized our understanding of 
the formation and evolution of galaxies by uncovering a population of high- redshift, dust- 
obscured systems that are forming stars at a tremendous rate (e.g. Small et al. 1997; Hughes 
et al. 1998; Blain et al. 2002). Submillimeter galaxies (SMG) are high- redshift analogues of 
the local ultraluminous infrared galaxies (ULIRGs) in terms of their CO, radio and infrared 
properties (Tacconi et al. 2008, 2010). The general importance of ultraluminous sources 
for galaxy evolution has been highlighted by recent Spitzer and Herschel results that show 
that thes e sources contribute s ignificantly to the total amount of star formation in the early 



universe ( iMagnelli et al. 20091 ). 



Until recently, observations at submillimeter wavelengths have been conducted mainly 
from the ground and/or from balloon experiments and covered only relatively small areas of 
the sky. The discovery rate of rare objects such as strongly lensed sources was therefore lim- 
ited. Nevertheless, due to selection effects, many of the known high-redshift submillimeter 
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sources ar e Kravitationally lensed, such as /■fi^6'F10214+47 24, which is amphfied by a factor 
of ~ 10 X (JRowan- Robinson et al. 199ll : lDownes et al. 19951 ). However, these hmited surveys 
found few strongly lensed sources. This situation changed dramatically with the 87deg^ 
deep field survey a t wavelengths of 1.4 and 2.0 mm taken with t he South Pole Teles cope 
( IVieira et al. 20101 ) and, since then, the Herschel space observatory (jPilbratt et al. 2010l ) has 
been able to map an enormously extended area of the sky at submillimeter wavelengths. In 
particular, the wid e-area Herschel Astrophysical Terah ertz Large Are a Survey or iJ-ATLAS 
JEales et al. 201oh with the SPIRE dOriffin et al. 201o[ ) and the PACS JPoglitsch et al. 201o[ ) 
instruments will ultimately map a total of about 570 square-degrees to around (or below) 
the confusion limit in five bands from 100 to 500 /im. In addition to uncovering hundreds 
of thousands of SMGs, the Herschel surveys will enable discovery of the brightest infrared 
sources in the universe, including many strongly len sed sources because of the high magni- 
fication bias in the submillimeter bands ( iBlain 19961 ). These surveys and samples of lensed 
sources are unique and will remain so for the foreseeable future. 

Negrello et al. (2007) predicted that most of the brightest 500- /im sources (i.e. those with 
S^oofim > 100 mJy) detected at high redshift by the Herschel space observatory will consist 
of strongly lensed submillimeter galaxies, a population of local star-forming galaxies - easily 
identified in shallow Sloan Digital Sky Survey (SDSS) images - a relatively small fract ion of 
radio-bright AGNs - easily identifiable in shallow radio surveys (JNegrello et al. 20101 ) - as 
well as nearby spirals. The ongoing surveys have confirmed this prediction with the detection 
of submillimeter-bright (S'soo/xm > 100 mJy) sources that have been shown, based on ground- 
based foll ow-up observations, to be gravitationally lensed SMGs in the distant (2; ~ 1.5 — 4) 
Universe ( jNegrello et al. 2010t IConley et al. 20111 ). The lensing boosts the sensitivity of 
observations and improves their spatial reso lution enabling de t ailed studies of the po pulations 
responsible for the bulk of the background (JBlain et al. 1999t iHopwood et al. 20101 ) . Follow- 
up observations with ground-based telescopes, especially by (sub)millimeter facilities are 
essential to further probe the physical characteristics of these outstanding sources. A recent 
and striking example of a lensed source d iscovered in the \)Te- H erschel era, in the submillimeter 
is the z = 2.3 galaxy SMM J2135-0102 JSwinbank et al. 201ol ) or the 'Cosmic Eyelash'. This 
typical SMG, which is amplified by a factor of 32.5 x, has been obser ved in exquisite detail 



enabling diagnosis of the properties of its molecular and atomic gas ( iDanielson et al. 2011 
Swinbank et al. 201 ll ). 



Here, we report ground-based follow-up observations of an exceptional source uncovered 
by the i/-ATLAS survey, HATLAS J142413.9+022304 (hereafter ID 141), including a mea- 
surement of its redshift and a derivation of the properties of its molecular and atomic gas 
content. 
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Throughout the pa per we use a ACDM cosmology with Hq = 71kms ^ Mpc ^, Qm 



0.27 and ^a = 1 - ^m flSpergel et ah 20031 ) 



2. Observations and results 

The source ID 141 is one of the brightest sources detected so far in the Herschel deep 
surveys. It is the strongest source peaking at 500 fim in the 100 deg^ currently surveyed 
in the i7- ATLAS program and yet found in the Herschel cosmological surveys. The source 
position as determined by Herschel is RA 14:24:13.9, Dec. +02:23:04 (J2000) and is located 
in the so-called 'GAMA 15h' field. Its flux densities, as measured with SPIRE, are in excess 
of 100 mJy in all bands with a spectral energy distribution (SED) that is rising from 250 
to 500 /im (Tabled]), indicating that the redshift of this exceptional source must be high 
{z>3). 

Imaging from the SDSS (DR7 - York et al. 2000) reveals a faint optical counterpart 
with the following magnitudes u = 23.45 ± 1.33, g = 22.99 ± 0.35, r = 22.06 ± 0.23, 
i = 20.79 ±0.12, and z = 20.35 ±0.40, that is located within 0.4" of the far-infrared position 
as measured by Herschel. The photometric redshift based on the optical magnitudes is 
z = 0.69 ± 0.13, which is clearly incompatible with the high redshift indicated by the far- 
infared photometry and mo re consistent with the hypothesis that this source is the lensing 



galaxy (JNegrello et al. 2010l ). 



2.1. 30-nieter Telescope 

ID 141 was subsequently observed at the IRAM 30-meter telescope on 2010 May 19 using 
the 117-channel bolometer array MAMBO-2 to measure the dust continuum flux density at 
1.2 mm. The observations were made using the on-off mode, wobbling by 32" in azimuth 
at a rate of 2 Hz. The target source was positioned on the most sensitive bolometer, and 
the correlated sky noise was determined from the other bolometers and subtracted from the 
'on-source bolometer'. The source was observed twice, each time for 10 min, yielding a final 
r.m.s. noise level of ~ 2 mJy. Standard sources (planets and secondary calibration sources) 
were used for absolute flux calibration. The data were reduce d using the standard scripts 



for on-off observation data in the MOPSIC software package (IZylka 19881 ). The measured 
flux density of ID 141 is 5*1.210111 = 36±2mJy, making it one of the high-redshift sources with 
the largest flux density measured to date at millimeter wavelengths. 

Together with the Herschel data, the 1.2-mm measurement indicates that the photo- 
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metric redshift is in the range 3.3 < 2;phot < 4.5 - following the method of fitting low-redshift 
templates to the far-infrared/millimeter available data described in Aretxaga et al. (2005) - 
a range containing 68% of the probability distribution that corresponds to a very fiat distri- 
bution due to the lack of photometric data between 500 /im and 1.2 mm, where the peak of 
the spectral energy distribution (SED) is lying. 



2.2. Plateau de Bure Interferometer (PdBI) 

Based on the redshift range estimated from the available photometric data of Herschel 
and MAMBO, we used the IRAM PdBI (with five antennas) to search for CO emission 
lines and determine the redshift of ID 141. This search was made possible by the recently in- 
stalled wide-band correlator WideX that provides a contiguous frequency coverage of 3.6 GHz 
bandwidth in dual polarization with a fixed channel spacing of 1.95 MHz. The observations 
were made in D configuration from 2010 May 23 to May 27 in a series of brief integrations 
with good atmospheric phase stability (seeing of 1") and reasonable transparency (PWV 
of 5-10 mm). Absolute fluxes was calibrated on MWC 349 when available or the quasar 
1502+1749 whose flux is regularly monitored at the IRAM PdBI. The accuracy of the flux 
calibration is estimated to be ~10% at 3 and 2 mm. 

We started the line search by tuning the central frequency of the band to 98.6 GHz 
(96.8-100.4 GHz) with the aim to cover in a series of observations the frequency range from 
104.0 to 86.0 GHz that is encompassing the redshift range (3.5 < z < 4.5) based on the 
probability distribution of Zphot- For this range in redshift, the expected CO transition is the 
4 — )■ 3. For each setting, the r.m.s. noise was la ~ 4mJy in 10 MHz channels after 20 min 
on-source. After five frequency set-ups, covering 18 GHz from 104 to 86 GHz, a strong and 
broad emission line was detected at 87.9 GHz with a peak intensity of ~ lOmJy and a zero 
intensity width of ~ 1000 kms"^ (Fig. [1] and Table [2]). 

To identify the transition of the CO emission line detected at 87.9 GHz and search for 
the next higher CO transition, we tuned the receivers to a central frequency of 108.5 GHz 
covering a redshift range that includes the most probable case where the 87.9 GHz line is 
the 4^3 transition of CO. A strong line was immediately apparent centered at 109.9 GHz, 
with a peak intensity of ~ 13 mJy and a broad line width comparable to that of the lower 
frequency line. The two lines at 87.9 and 109.9 GHz are identified with the 4-^-3 and 5 — )■ 4 
transitions of CO (Figs. [T]&|2]). These identifications yield a redshift oi z = 4.243 ± 0.001 
for ID 141. This is the first 'blind' redshift determination using the PdBI. 

In addition to the CO emission lines, the 3-mm dust continuum emission was detected 
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at a level of 1.6 ± 0.1 mJy, averaged over the entire line-free range covered at 3 mm, i.e. over 
~ 19 GHz. Tabled] lists the continuum flux densities derived from this 3 mm spectral survey 
for three different wavelengths (2.75, 3.00 and 3.29 mm) each averaged over ~ 4 GHz. 

Subsequent observations were performed on June 29 and September 11 2010, using the 
compact D configuration under moderate weather conditions (PWV of 5-10 mm, 100 < 
Tgys/K < 300) to search simultaneously for the emission lines of C0(7 — t- 6) and [Cl](^P2 —^ 
Pi) redshifted to 153.876 and 154.339 GHz. Both lines were detected with relatively good 
signal-to-noise ratio together with the 2-mm dust continuum at a level of 9.7 ± 0.9 mJy 
(Fig. [2]). The detection of the [Cl](^P2 —^ Pi) emission line prompted us to check for the 
presence of the [Cl](^Pi — ^ Pq) transition that is redshifted to 93.87 GHz and is included in 
the frequency sweep done from 104.0 to 86.0 GHz. Due to the low signal-to-noise ratio of 
the sweep, the [Cl]('^Pi — '^ Pq) line is detected with a poor signal-to-noise ratio and a line 
intensity of 2.8 ± 0.9 Jykms^^. 

The parameters of the CO and [Cl] lines observed towards ID 141 are listed in Table [2l 
Figure El summarizes the results based on the PdBI observations in a series of maps displaying 
both the continuum and line emission. The source, which is clearly not multiple, is barely 
resolved in the PdBI data at the best available beam of 3.6" x 3.0" at 2 mm. In the 2-mm 
continuum map, the integrated flux density of 9.7 ± 0.9 mJy is higher than the emission 
peak of 5.5 mJy, and a u-v fit yields a deconvolved size of 3.9 ± 0.4 x 2.4 ± 0.4arcsec^ (or 
(27 X 16)//iLkpc2) at an angle of -25 ± lOdeg (NE). 



2.3. SMA 

ID 141 was observed in the continuum at 880 fira with the Submillimeter Array (SMA) 
in its compact configuration with 8 antennas on 2010 June 15 under moderate weather 
conditions. The observations on source amounted to 2.4 hours. The blazar 3C279 was 
utilized as a bandpass calibrator, and Titan as the absolute flux calibrator. The resulting 
map is shown in Fig. [31 The source, which is well detected in the continuum, appears to 
be relatively small with respect to the PdBI source size and is resolved at the 2" resolution. 
The map r.m.s. noise is 1.75 mJy, with a peak of ~ 50 mJy. The integrated flux density 
is 90 ± 2 mJy. To within the limits of the observations, ID 141 appears to be elongated 
in the North- West direction, in agreement with the 2-mm measurements of the PdBI. The 
deconvolved size at 0.88 mm is 2.4 ± 0.1 x 1.3 ± O.larcsec^ (or (16 x 9) / fii^kpc^) at an 
angle of —32 ± 4deg (NE). The peak position of ID 141, as derived from the 0.88-mm map, 
corresponds to RA 14:24:13.98, Dec. +02:23:03.45 (J2000), which is slightly offset by 1.3" 
from the position determined by Herschel (see Fig. [3]). Clearly higher angular resolution 



observations (both in the sub/miUimeter and the optical) are required to further explore the 
structure of ID 141 and its relation with the lensing galaxy. 



2.4. APEX 

Observations of the redshifted 1900.538 GHz C [II] line emission in ID 141 were made 
with the APEX 12-meter telescope in 2010 July in a series of two observing sessions, using 
the refurbished dual-color FL ASH-!- MPIfR PI r eceiver. Its 345-GHz channel now operates 



wideband IRAM 2SB mixers (JMaier et al. 20051 ). The observations took place on July 7 



(where the source was observed for 86 minutes on-time) and July 25 (for 54 minutes on- 
time) during excellent weather conditions with zenith precipitable water vapour of 0.7 and 
0.4 mm. The spectra were taken in double-beam chopping mode, with a wobbler throw of 
50" at 1.5 Hz. Pointing was established on nearby Mars and Saturn. At the frequency of 
362.45 GHz, the beamsize of APEX is 17.8". 

The 4-GHz wide IF of the mixer was connected to a set of fast Fourier-transform spec- 
trometers. In the first observing session three spectrometers from the CHAMP -l-FFT back- 
end array, each 1.5 GHz wide, were combined to process, with overlap, a total of 3.4 GHz. 
In the second run - during their commissioning - two of the new 2.5-GHz-wide XFFT spec- 
trometers were connected to cover the full 4 GHz wide IF band (with 500 MHz of overlap), 
thereby providing 3300 km s~^ instantaneous velocity coverage. 

The [C II] emission line in ID 141 was readily detected in both observations. For data 
processing only linear offsets were removed. The combined spectrum has an r.m.s. noise 
of 1 mK. The line is detected with good signal-to- noise with a peak flux density of 4 mK 
in T* and a line shape that is comparable to that of the CO emission lines (Fig. [2]). The 
total integrated intensity of the [Cll] line is 2.64 ± 0.41 Kkms~^ or 107 ± 17 Jykms~^ (Ta- 
ble ED - at 362 GH z, the gain of the APEX is 40.7 Jy/K in converting T* to flux density 



(JGiisten et al. 20061 ). Because of the excellent weather conditions, the absolute calibration 
uncertainties are expected to be of order 10-20%. 

The continuum at 850 /xm was measured using the LABOCA 870-/im bolometer camera 
on 2010 September 20 for a total of about 10 minutes under excellent weather conditions. 
Planets as well as secondary calibrators were used for absolute flux calibration. The source 
was detected with a flux density of 102.1 ± 8.8 mJy, in good agreement with the SMA 
measurements (Table [1]). It is noted that the [ClI] contribution to the 60 GHz LABOCA 
passband amounts to only 2.1 mJy or 2%, which is smaller than the uncertainties of the 
current 850-/im continuum measurement. 
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Discussion 



We have detected and/or mapped at high significance three high-J emission hnes of 
CO and the fine-structure emission hnes of atomic ([Cl]) and singly ionized ([ClI]) carbon 
towards ID 141. Together with the detection and mapping of the far-infrared dust continuum 
emission, these results are used to constrain the properties of the molecular and atomic gas 
in this high-redshift galaxy. The detection of the CO emission lines firmly establishes the 
redshift of ID 141 at z = 4.243 ±0.001, a value within the redshift range derived from the far- 
infrared/millimeter photometric measurements. This redshift corresponds to a luminosity 
distance Dl = 3.9 x lO'^Mpc and a scale of 6.9 kpc/". 



3.1. Far-infrared spectral energy distribution 

Figure H] displays the SED of ID 141 - based on the flux densities listed in Table [1] 
- where it is compared with the SED of M82. Both SEDs are comparable from the rest- 
frame sub-millimeter to far-infrared wavelengths, both peaking at a rest wavelength of ~ 
90 /im. We fit the photometric data points of ID 141 with a single-temperature modified 
blackbody spectrum (using the optically thin approximation) by performing a x^ fit on 
both the dust temperature (Tdust) and the dust emissivity index (/3) as outlined in Beelen 
et al. (2006). We derive (3 = 1.7 ± 0.1, Tdust = 38 ± 1 K and a far-infrared luminosity of 
Lfir ~ 8.5 ± 0.3 X 10^^ /^L ^ L0, where /xl is the lensing magnification factor. To estimate the 
dust mass, we adopt for the mass absorption coefficient k, a value of 0.4cm^g~^ at 1200 /xm 
(e.g. Beelen et al. 2006). From the above dust temperature and far- infrared luminosity, we 
derive Mdust ~ 8.9 x IO^^Ul^Mq. Since the dust emission in luminous infrared starbursts 
can be optically thick around 100 /zm, we also modelled the SED of ID 141 following the 
analysis described in WeiB et al. (2007). A single component optically thick dust model 
(adopting f3 = 1.7) fits the photometric data equally well as the optical thin approximation, 
although with a higher dust temperature (Tdust ~ 58 K) and a lower dust mass (Mdust ~ 
5.5 X lO^yUL^ M0). The resulting fit is shown in Figure HI In this model, the dust becomes 
opaque (r = 1) at a rest-wavelength of 150 yum (corresponding to an observing frequency 
of 380 GHz) and the appar ent equivalent ra dius - the true galaxy radius for an unlensed. 



face-on, filled circular disk fjWeiB et al. 20071 ) - is ro ~ 720 fi^ ^ pc 



The global star-formation rate (SFR) can be estimated from the far-infrared luminosity 
using the relation SFR ^ 6mf (Tfir/10^° L fn) M^ yr~^ where (5m f is a function of the present 



mass composition of the stellar population (jOmont et al. 200ll ). For a range in metallicities 



initial mass functions and starburst ages of 10 to 100 Myr, 6mf ranges in between 0.8 and 2. 
Under the assumption that the far-infrared luminosity of ID 141 is dominated by starburst 
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activity, we estimate 0.7 < SFR < 1.7 x lO^/i^^ MqJi^^, which is an upper hmit in case 
there is an AGN contribution to the far-infrared luminosity. 



Using the FIRST radio survey (JBecker et ah 19951 ). the 1.4-GHz flux density of ID 141 



is measured to be S'i.4ghz = 570±160 /iJy, which shows that ID 141 is not a radio-loud source. 

For ID 141, the far-infrared/radio parameter q = logio [(-Fir/3.75 x 10^^ W m~^)/(S'i.4GHz/W m"^ Hz~^)] 



(ICondon 19921 ) is estimated to be g ~ 2.1 ±0.1, after i^-correcting and assuming a canonical 
radio slope for a star-forming galaxy a = —0.75, where 5*1^ oc z/°. This value is compati- 
ble with the median value g ~ 2.4 derived for starbursts samples selected from BLAST or 
the Herschd GOODS-North surveys (Ivison et al. 2010a, 2010b; Jarvis et al. 2010). This 
suggests that star-formation activity is reponsible for most of the far-infrared luminosity of 
ID 141. 

Typical amplification factors of 10 — 30 have been estimated for source s similar to ID 141 



with 500 /im flux density above 100 mJy uncovered in the i7- ATLAS survey (JNegrello et al. 20101 ). 
We therefore adopt in this paper a magnification range of 10 < /xl < 30 for ID 141. Under 
this assumption, the far-infrared luminosity, dust mass and global star formation of ID 141 
are estimated to be in the range: 2.8 to 8.5 x 10^^ L© for Lpm, 3 to 9 x 10^ Mq for M^ust, and 
200 to 1700 Mq yr~^ for SFR. These values are typical for the population of dusty , gas-rich 



and luminous SMGs found at high redshift (JGreve et al. 20051 : iTacconi et al. 2010l ) 



3.2. The CO emission lines 

The CO lines of ID 141 are broad (AVfwhm ~ 800kms~^ - see Table [2]) and display a 
clear, asymmetrical and double-peaked profile where the blue part of the line is stronger than 
its red counterpart (Fig. [2]) . A similar profile is seen in the [C II] emission line. The blue- and 
red-shifted peaks of the lines are found from Gaussian fits at velocity offsets of ±220kms~^, 
respectively, each with a FWHM of ~ 400kms~^. The broad line width of ID 141 com- 
pares with the median FWHM (780 ±320 km s~^) of high-redshift SMGs and double-peaked 
asymmetrical line s hapes have also been reported for other luminous submillimeter galaxies 



(jGreve et al. 20051 : ISwinbank et al. 20101 ). Such double-peaked line profiles indicate orbital 
motions under the influence of gravity either in the form of a disc or a merger. However, 
the current angular resolution of the observations presented in this paper does not allow to 
distinguish between the two scenarios and attempts to find positional offsets between the 
blue- and red-shifted peaks remained inconclusive. Higher resolution observations in either 
CO or [Cll] are required to further examine the structure of the velocity gradient in ID 141. 

From the line parameters listed in Table [2] and following the definition for line luminosity 
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given by Eq. (3) in Solomon et al. (1997), the CO line luminosities are L' = 3.34 ± 0.40, 
3.71 ±0.45 and 0.94±0.20 x 10" fi^^ K kms"^ pc^ for the 4 ^ 3, 5 ^ 4 and 7 ^ 6 transitions, 
respectively. This corresponds to CO(J = 5 ^ 4)/C0(J = 4—^3) and CO(J = 7 — )■ 
6)/C0(J = 5 — > 4) line brightness temperature ratios of r54 = 1.1 ±0.1 and rys = 0.24±0.04, 
indicating that the higher CO transitions are not thermalized (see also values of Lqq in 
Table ED. 

To describe the properties of the molecular gas in more detail, we have analyzed 
the CO data usi ng a spherical, single-component, large velocity gradient (LVG) model 
(IWeifi et al. 20071 ). As in WeiB et al. (2007), we adopt a CO abundance per velocity gradient 
of [CO]/{dv/dr) = 10~^pc (kms^^)~^ that best fits the high- J lines and implies low enough 
opacities to reproduce the fluxes of the two lower CO transitions. The results of the best-fit 
model, which are displayed in Fig. correspond to an H2 density of nnj = 10^'^ cm~^, a gas 
temperature of Tkin = 40 K and an equivalent (magnified) radius of tq = 680 pc. These val- 
ues are comparable to those derived for other hi gh-redshift SMGs such as GN20 aX z = 4.04 
JCariUi et al. 201o[ l and SMM J16359 aX z = 2.5 JWeifi et al. 2005bl ). as shown in Fig. El Typ- 
ical degeneracies between the temperature and density in such LVG models are discussed in 
WeiB et al. (2007). Additional CO transitions, in particular the low- J transitions, should be 
observed to evaluate more precisely the density and the temper ature of the molecu lar gas and 
check whether an excess of cold molecular gas could be present (jivison et al. 201 if ). Based on 
this model, the expected C0(1— 0) line luminosity is -^co(i-o) ~ 4.3x10^^/Xl^ Kkms~^ pc^ or 



^CO(l-O) 



2.1 X 10 /xlL©. Assuming a conversion factor of aco = 0.8 Mq (Kkms pc 



from L^Q^QN to Mh,, appropriate for starburst galaxies as derived from observations of 
local ULIRGs (JDownes fc Solomon 19981 : iTacconi et al. 20081 ) . yields a molecular gas mass 
Mhj = 3.5 X 10^^ /iL^ Mq. After correction for the amplification, the above values are com- 
parable to the median gas mass (3 x 10 ^° M©) derived for luminous SMGs using the same 
conversion factor aco (jGreve et al. 20051 ) and the gas masses of (1 — 3) x 10^° M© measured in 
the i/-ATLAS high -redshift lensed SM Gs, SDP.81 and SDP.130, discovered by the Herschel 
space observatory (JFrayer et al. 20101 ). 



3.3. The Carbon fine-structure atomic lines 



3.3.1. The [C II] emission line 

From the total integrated intensity, the [ClI] line luminosity for ID 141 is Lcii = 6.1 ± 
0.9 X IO^'^/Xl^Lq. Correcting for the amplification factor, the value of Lqh for ID 141 is 
therefore comparable to that of other high-redshift sources, such as: SP SS J1148±5251 
{z = 6.42) with 4.4 ± 0.5 x 10^ Lq flMaiolino et al. 20051 : IWalter et al. 2009h . BR 0952-0115 
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4.43) with 4.6 ± 0.7 x lO^L^ flMaiolino et al. 2009h . BRI 1335-0417 



4.40) with 



16.4 ±2.6 X IO^Lq JWagg et ah 201oh. SMM J2135-0102 (the 'Cosmic Eyelash' at z = 2.3) 



with 5.5 ± 1.3 X W^Lq fllvison et ah 2010cr). or the submiU imeter galaxy, LESS J033229.4 



{z = 4.76) with 1.02 ±0.15 x 10^° Lq flDe Breuck et al. 20111 ). The [ClI] line luminosities for 
these high-redshift infrared-luminous galaxies are, however, lower than the most of the values 
derived for a sample of 14 galaxies at redshifts 1 < z < 2 that have bright [Cll] emission lines 
yielding luminosit ies in the range 8 x 10^ < Lq^JLq < 10^^ (JHailey-Dunsheath et al. 2010l : 
Stacev et al. 201oh . 



For ID 141, we derive a ratio L[Cii]/-^fir = (7.3 ± 1.31) x 10 "^. This is higher than for 
most of the high-redshift sources with similar infrared l uminosity and reds hift but comparable 
to BRI 1335-0417 where the ratio is (5.3±0.8) x 10"^ JWagg et ah 201oh . ID 141 follows the 
trend reported for high-redshift luminous galaxies in the relation of L[Cii]/-^fir as a function 
of the far-infrared luminosity Lpm, with a ratio of about one order of magnitude lower than 
for normal galaxies but comparable to ULIRGs values (Fig. |6]). 

Recently Gracia-Carpio et al. (2011) showed that the ratio between the far-infrared 
luminosity and the molecular gas mass, Lfir/Mh2, appears to be a better choice from a 
physical point of view to investigate the relative brightness of the far-infrared fine-structure 
lines in galaxies. This ratio, which removes the uncertainty of the amplification factor, also 
appears to better reflect the properties of the gas than Lfir alone. Gracia-Carpio et al. 
(2011) found that there is a threshold at Lfir/Mh2 ~ SOLqMq^ above which galaxies tend 
to have weaker [C II] emission with respect to their far-infrared luminosity. This threshold is 
similar to the one separating the star-formati on relations between non- or weakly-in teracting 
star- forming galaxies and luminous mergers (JGenzel et al. 2010t iDaddi et al. 20101 ). With a 
ratio Lfir/Mh2 ~ 23OL0M0\ ID 141 follows the trend in the L[Cji]/I/fir vs. Lfir/Mhj 
plane reported in Gracia-Carpio et al. (2011), indicating that its ionization parameter must 
be high (~ 10~^) and that ID 141 could be a luminous merger. 

Another sensitive probe of the physical conditions within the interstellar medium is the 
ratio Lcii/i^co{i-o)- For ID 141, we find from the measured [ClI] and the estimated CO(l-O) 
luminosities, Lcii/i^co(i-o) ~ 2900. Using the Photon-Dominated Regioi i (PDR) model and 
diagnostic diagram comparing the ratios Lc„/-^fir versus Lco(i^o)/-^fir flKaufman et al. 19991 ) 
indicates that, for ID 141, the gas is dense with n ~ 10^ cm~'^ and that the far-UV radiation 
field illuminating the PDR, Go, is a few 10^ times more intense than that in our Galaxy. Both 
these values po int to physical con ditions comparable to those found in the central regions of 
local ULIRGs (Istacev et al. 201ol ). 
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3.3.2. The [Cl] emission line 

The [CI](=^P2 -^ Pi) line flux implies Lci = 0.8 ± 0.3 x 10 Vl ^ L© (Table E]). The [Cl] 
luminosity in ID 141, once corrected for the amplification, is therefore lower than the [Cl] 
luminosities de rived from the (^Py, —'^ Pi ) transitio n for other high-reds hift sources such as 
the Cloverleaf JWeifi et al. 20031 : IWeifi et al. 2005al) or PSS 2322+1944 JPetv et al. 2004} ) or 
other high-redshift galaxies or quasars (jWalter et al. 201 ll ). Assuming an excitation tem- 
perature (Tex) equal to the temperature of the dust and gas, i.e. Tex ~ 40 K, the mass 
of neutral carbon in ID 141 ar nounts to Mct = = 6.5 ± 1.1 x 10"^ h^^Mq, a result that is 
only weakly dependent on Tex (jPety et al. 2004f l. Compared to the mass of molecular gas, 
the derived mass of [Cl] implies a neutral carbon abundance relative to H2 of [CI]/[H2] ~ 
1.8 X 10~^. The similarity of this abundance w ith the value of 2.2 x 10~^ found for Galac- 
tic dense molecular clouds (JFrerking et al. 19891 ) is indicative of a substantial enrichment in 
ID 141. In addition, the LcTfi-m/TpTR lumino sity ratio is also a sensitive measurement of Go 
kaufman et al. 1999t ICerin fc Phillips 200ol ). For ID 141 Lmn-i^^/Lym ~ 5 x 10"^ close to 



the luminosity ratio found in other high-redshift sources ( jWalter et al. 201 if ) and indicating 
a high UV field strengh of Go ~ 3 x 10^, compatible with the analysis based on the [ClI] 
emission line. 

The measurements available for ID 141 allow an estimate of the relative cooling of the 
carbon and CO emission lines. The upper transition of [Cl] is stron ger than the lowe r 
transition by a factor of ~ 2 (Table [2]), as in the case of the Cloverleaf (jWeiB et al. 2005a[ ). 
The total luminosity of ID 141 in [Cl] is therefore ^ 1.2 x IO^/Il^L©. Based on the LVG 
model described above, the total CO line luminosity in ID 141 is estimated to be Leo ~ 
8.4 X 10^ /iL^ Lq. The CO/[Cl] luminosity ratio is thus ~ 7.0 and the cooling in the CO and 
[C I] lines represents 10~^ of the far-infrared continuum, a factor of 7 lower than the main 
cooling line of the gas, i.e the [C II] transition. 



Conclusions 



We have used four ground-based (sub) millimeter facilities - the IRAM 30- meter and 
PdBI, the SMA and APEX - to follow up the exceptional source ID 141 that is the strongest 
500-/im riser yet discovered in the Herschel cosmological surveys. The detection of three 
bright CO emission lines using the PdBI confirms that the source is at high redshift {z = 
4.243). The strength of the continuum and the intensities of the molecular and atomic 
emission lines of ID 141 suggest the source is gravitationally lensed, in line with what was 
already found for other sub- millimeter brigh t , high-redshift, galax ies detected with the Her- 
schel space observatory (JNegrello et al. 20071 : iNegrello et al. 2010l ). The magnification factor 



-14- 



is still unknown but is likely to be in the range 10 < /xl < 30 as found for other gravita- 
tionally lensed SMGs uncovered in the i7-ATLAS survey. The properties of the molecular 
and atomic gas in ID 141 indicate properties of a luminous (LpiR ~ 3 — 8 x 10^^ -^o), dense 
(n ~ lO^cm"^) and warm (Tkin ~ 40 K) starburst galaxy that are comparable to those de- 
rived for local starburst galaxies or high-redshift submillimeter galaxies. The asymmetric 
double profiles of the carbon and CO emission lines are indicative of orbital motions in a 
disc or a merger. The source is barely resolved in the present data indicating a physical size 
that is smaller than 14/yUL kpc, corresponding to an upper limit of 2". 

The measurements of the carbon lines ([Cl] and [ClI]) in ID 141 enable a first study 
of the gas associated with the PDR in ID 141 and analysis of the relative cooling of the 
carbon and CO emission in this galaxy, that is dominated by that of C^. The [ClI]-to- 
FIR luminosity ratio, (7.3 ± 1.3) x 10"'^, is larger than the one found for other galaxies of 
comparable luminosity and redshift, and consistent with recent findings that the scatter in 
the Lcu/Lpi^ versus Lpm diagram is indeed large. With a ratio Lfir/Mhj ~ 23OL0M0\ 
ID 141 follows the trend reported by Gracia-Carpio et al. (2011) in the Lcii/i^FiR versus 
Lfir/Mmh, plane, suggesting both a high ionization parameter and the fact that ID 141 
could be a merger. 

The present observations underline the importance of observations done at frequencies in 
between 85 and 360 GHz to study the properties of galaxies at redshift z > 4. In particular, 
the lower frequencies (in the 3 and 2 mm windows) are essential to probe the peak of the CO 
emission that, for such starburst galaxies, lies around Juppcr 4 to 5, as well as the two emission 
lines of [C I] or molecules other than CO such as water. Follow-up observations of very-high- 
redshift, obscured galaxies found in future deep surveys, especially the less luminous objects, 
will need sensitive interferometers operating at frequencies below 350 GHz such as ALMA, 
the EVLA or the future upgrade of the Plateau de Bure interferometer (NOEMA) in order 
to fully explore the properties of the molecular and atomic gas in starburst galaxies that 
were active when the universe was about 1 Gyr old. 

The results described in this paper are based on observations obtained with Herschel, an 
ESA space observatory with science instruments provided by European-led Principal Investi- 
gator consortia and with important participation from NASA. US participants in i/- ATLAS 
acknowledge support from NASA through a contract from JPL. The ground-based follow-up 
observations were obtained at the following facilities. The 30-meter telescope and the PdBI 
of IRAM that is funded by the Centre National de la Recherche Scientifique (France), the 
Max-Planck Gesellschaft (Germany) and the Instituto Geografico Nacional (Spain). APEX is 
a collaboration between the Max-Planck-Institut fiir Radioastronomie, the European South- 
ern Observatory, and the Onsala Space Observatory. The Submillimeter Array (SMA) is a 
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joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. Part of this study was supported by financial contribution from the 
agreement ASI-INAF 1/009/10/0 and by CONACyT grant 39953-F. 

Facilities: Herschel (SPIRE), IRAM 30-meter, IRAM PdBI, SMA, APEX. 



REFERENCES 

Aretxaga, I., Hughes, D.H., Dunlop, J.S. 2005, MNRAS, 358, 1240 

Becker, R.H., White, R.L., Helfand, D.J. 1995, ApJ, 450, 559 

Beelen, A., Cox, P., Benford, D. et al. 2006, ApJ, 642, 694 

Blain, A.W. 1996, MNRAS, 293, 1340 

Blain, A.W., Kneib, J.-R, Ivison, R.J., Smail, I. 1999, ApJL, 512, L87 

Blain, A.W., Smail, I., Ivison, R. J., Kneib, J.-P., Prayer, D. T. 2002, Physics Reports, 369, 
111 

Carilh, C.L., Daddi, E., Riechers, D. et al. 2010, ApJ, 714, 1407 

Carral, P., Hollebach, D.J., Lord, S.D. et al. 1994, ApJ, 423, 223 

Colbert, J.W., Malkan, M.A., Clegg, RE. et al. 1999, ApJ, 511, 721 

Condon, J. J., ARAA, 30, 572 

Conley, A., Cooray, A., Vieira, J.D. et al. 2011, ApJ 732, L35 

Daddi, E., Elbaz, D., Walter, F. et al. 2010, ApJ, 714, LI 18 

Danielson, A.L.R. et al. 2011, MNRAS, 410, 1687 

De Breuck, C, Maiolino, R., Caselli, P. et al. 2011, A&A, 530, L8 

Downes, D., Solomon, P.M., Radford, S.J.E. 1995, ApJ, 453, L65 

Downes, D. & Solomon, P.M. 1998, ApJ, 507, 615 

Eales, S., Dunne, L., Clements, D. et al. 2010, PASP, 122, 499 



-16- 

Frayer, D.T., Harris, A.L, Baker, A.J. et al. 2010, ApJ, 726, 22 

Frerking, M.A., Keene, J., Blake, G.A., & Phillips, T.G. 1989, ApJ, 344, 311 

Genzel, R., Tacconi, L.J., Gracia-Carpio, J. et al. 2010, MNRAS, 407, 2091 

Gerin, M. & Phillips, T.G. 2000, ApJ, 537, 644 

Gracia-Carpio J., Sturm E., Hailey-Dunsheath S. et al. 2011, ApJ, 728, L7 

Greve, T.R., Bertoldi, F., Small, I. et al. 2005, MNRAS, 359, 1165 

Griffin, M.J., Abergel, A., Abreu, A. et al. 2010, A&A, 518, L3 

Glisten, R., Nyman, L.A., Schilke, P. et al. 2006, A&A, 454, L13 

Hailey-Dunsheath, S., Nikola, T., Stacey, G.J., Oberst, T.E. Parshley, S.C., Benford, D.J., 
Staguhn, J.G., Tucker, C.E. 2010, ApJ, 714, L162 

Hopwood, R. et al. 2010, ApJL, 716, L45 

Hughes, D.H., Serjeant, S., Dunlop, J. et al. 1998, Nature, 394, 241 

Ivison, R.J., Alexander, D.M., Biggs, A.D. et al. 2010a, MNRAS, 402, 245 

Ivison, R.J., MagneUi, B., Ibar, E. et al. 2010b, A&A, 518, L31 

Ivison, R.J., Swinbank, A.M., Swinyard, B. et al. 2010c, A&A, 518, L35 

Ivison, R.J., Papadopoulos, P.P., Small, I. et al. 2011, MNRAS, 411, 1913 

Jarvis, M., Smith, D.J.B., Bonfield, D.G. et al. 2010, MNRAS, 409, 92 

Kaufman, M.J., Woffire, M.G., Hollenbach, D.J., & Luhman, M.L. 1999, ApJ, 527, 795 

Luhman, M.L., Satyapal, S., Fischer, J. et al. 1998, ApJ, 504, Lll 

Luhman, M.L., Satyapal, S., Fischer, J. et al. 2003, ApJ, 594, 758 

Maier, D., Barbier, A., Lazareff, B. & Schuster, K.F. 2005, Proceedings of the 16**^ Interna- 
tional Symposium on Space Terahertz Technology, pp. 428-431 

Maiolino, R., Cox, P., Caselh, P. et al. 2005, A&A, 440, L51 

Maiolino, R., Caselh, P., Nagao, T., Walmsley, M., De Breuck, C. & Meneghetti, M. 2009, 
A&A, 500, LI 



■-; ^^^; 



-17- 

Magnelli, B. et al. 2009, A&A, 496, 57 

Malhotra, S., Kaufman, M.J., HoUenbach, D. et al. 2001, ApJ, 561, 766 

Negrello, M., Perrotta, F., Gonzalez-Nuevo, J. et al. 2007, MNRAS, 377, 1557 

Negrello, M., Hopwood, R., De Zotti, G. et al. 2010, Science, 330, 800 

Omont, A., Cox, R, Bertoldi, F. et al. 2001, A&A, 374, 371 

Pety, J., Beelen, A., Cox, P. et al. 2004, A&A, 428, L21 

Pilbratt, G.L. et al. 2010, A&A, 518, LI 

Poglitsch, A., Waelkens, C, Geis, N. et al. 2010, A&A, 518, L2 

Riechers, D.A., Capak, P.L., Carilli, C.L. et al. 2010, ApJ, 720, L131 

Rowan-Robinson, M., Broadhurst, T., Oliver, S.J. et al. 1991, Nature, 351, 719 

Small, I., Ivison, R.J., Blain, A.W. 1997, ApJL, 490, L5 

Solomon, P.M., Downes, D., Radford, S.J.E., Barrett, J.W. 1997, ApJ, 478, 144 

Spergel, D. N., et al. 2003, APJS, 148, 175 

Spinoglio, L., Malkan, M.A., Smith, H., Gonzalez- Alfonso, E., Fischer, J. 2005, ApJ, 623, 
123 

Stacey, G.J., Hailey-Dunsheath, S., Ferkinhoff, C. et al. 2010, ApJ, in press 

Swinbank, A.M. et al. 2010, Nature, 464, 733 

Swinbank, M. et al. 2011, submitted to ApJ 

Tacconi, L.J., et al. 2008, ApJ, 680, 246 

Tacconi, L.J., et al. 2010, Nature, 463, 781 

Unger, S.J., Clegg, P.E., Stacey, G.J. et al. 2000, ApJ, 355, 885 

Valtchanov, I., Virdee J., Ivison, R.J. et al. 2011, MNRAS, in press flarXiv:1105.3924l) 

Vieira, J.D., Crawford, T.M., Switzer, E.R., et al. 2010, ApJ, 719, 763 

Wagg, J., Carilh, C.L., Wilner, D.J. et al. 2010, A&A, 519, LI 



-18- 

Walter, F., Riechers, D., Cox, P. et al. 2009, Nature, 457, 699 

Walter, F., WeiB, A., Downes, D., Decarli, R., & Henkel, C. 2011, ApJ, 730, 1^ 

WeiB, A., Henkel, C, Downes, D., & Walter, F. 2003, A&A, 409, L41 

WeiB, A., Downes, D., Henkel, C, & Walter, F. 2005a, A&A, 429, L25 

WeiB, A., Downes, D., Walter, F. & Henkel, C. 2005b, A&A, 440, L45 

WeiB, A., Downes, D., Neri, R. et al. 2007, A&A, 467, 955 

York, D.G. et al. 2000, AJ, 120, 1579 

Zylka, R. 1988, MOPSI Users Manual, (IRAM: Grenoble) 



This preprint was prepared with the AAS E^TJtjX macros v5.2. 



-19- 



Table 1: Photometric Measurements of ID 141 



Wavelength 


Flux Density 


Facility 


fira 




mJy 




250 




115±19 


SPIRE^ 


350 




192±30 


SPIRE^ 


500 




204±32 


SPIRE^ 


870 




102±8.8 


LAB0CA2 


880 




90±5 


SMA 


1200 




36±2 


MAMBO-23 


1950 




9.7±0.9 


PdBI 


2750 




1.8±0.3 


PdBI 


3000 




1.6±0.2 


PdBI 


3290 




1.2±0.1 


PdBI 


Note - ^ 


Herschet ^ APEX: 


3 30-meter. 



Table 2: Parameters of the CO, [Cl] and [ClI] emission lines observed towards ID 141 



Line 



z/obs Peak Int. AVfwhm I 

[GHz] [mJy] [kms-i] [Jy km s^^] 



L L' 

[IO^Lq] [10i°Kkms-ipc2] 



CO (4-3) 


87.929 


9.2±0.7 


760±70 


7.5±0.9 


l.OiO.l 


33.4±4.0 


CO (5-4) 


109.911 


13.6±1.1 


890±80 


13.0±1.6 


2.3±0.3 


37.1±4.5 


CO (7-6) 


153.876 


7.7±1.2 


790±130 


6.5±1.4 


1.6±0.3 


9.4±2.0 


[CI](3Pi-3Po) 


93.870 


3.0±1.0 


790t 


2.8±0.9 


0.42±0.15 


10.9±3.5 


;ci](3p2-3Pi) 


154.339 


4.0±1.0 


790±130 


3.4±1.1 


0.8±0.3 


4.9±0.9 


[CII](2P3/2-'Pl/2) 


362.45 


147±20 


690±80 


107±17 


61.6±9.8 


28.1±4.46 



Note - The line luminosities are not corrected for lensing magnification. The amplification 
factor is supposed to be in the range 10 < /iL < 30 - see discussion in Sect. 3. ^ Fixed 

linewidth. 
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Fig. 1. — The full frequency sweep done at the PdBI to search for CO emission lines in ID 141, 
showing the detections of CO(J = 4 — )■ 3) and CO(J = 5 — > 4). The 4 — )■ 3 transition was 
detected first after covering 18 GHz in 5 frequency set-ups. The 5 ^^ 4 emission line was 
detected thereafter in one setting, confirming the redshift of the source at z = 4.243. The 
line shows the continuum at a wavelength of 3 mm that is clearly detected and is increasing 
from 1.2 ± 0.1 to 1.8 ± 0.3 mJy from the lowest to the highest frequencies - as determined 
from the averaged line- free channels - see Tabled]). The position of the tentatively detected 
[Cl]('^Pi — '^ Pq) emission line is indicated. 



-21 - 




_ 


20 


) : 


15 


- 


10 


1 - 


5 


'V 





' 





C0(7-6) - 200 




- 100 - 



100 




-2000 



1000 







1000 2000 



-2000 



1000 



1000 2000 



Velocity Offset (km s ) 



3) and CO(J 

2 



5 — )> 4) (upper panels), the 



Fig. 2.— Spectra of the CO(J = 4 ^ 

CO(J = 7^6), [CI](=^P2 -^ ^i) and [CU](^P3/2 -^ P1/2) emission hnes (lower panels) 
towards ID 141. The continuum was subtracted for all spectra. The velocity scale corre- 
sponds to the frequencies (z/obs) listed in Table |2] for the CO and [ClI] emission lines (i.e. for 
the redshift of 2; = 4.243). Gaussian fits to the CO(J = 7 — )► 6) and [Cl] emission lines are 
shown as dashed curves whereas the full line shows the composite fit. 



-22- 



jDntiXjuum a: 3.0mrn (PdBl) 



o. 



plntlnufCHon at 2.0 mm (PdBl}..-- 



UptOQClt 



(Conti/iuum at O.SBmmg (;5MA) . 
^-0 oo\ 



O Q ^'^ 



-10 




o C 



-- O °Cc 





^ 







B£Atvi - " 



o 

O -, 



• 
BEAM - 



C0(J=A^3) (PdBi; 



^6 



C0(J = 5-4) (PdBi; 



■■€0(J = 77-6))& CI(¥2-^P,)(4b5bI) 



10 - 



'd' 



o 






o 



o 



L 




9 



10 -10 

RA - offset (") 



10 -10 

RA - offset (") 



10 -10 

RA - offset (") 



Fig. 3. — Upper panels Maps of the dust continuum emission at 3, 2 and 0.88 mm towards 
ID 141. Lower panels: Velocity-integrated maps of the three CO and [Cl](^P2 — ^ -Pi) emis- 
sion hues towards ID 141. All the data were obtained with the PdBl except the 0.88-mm 
continuum emission that was measured with the SMA. The contours are shown in steps of 
2o", starting at ±2o" and the synthesized beams are shown in the right lower corner of each 
panel. The la values are: 0.09, 0.3 and 1.75 mJy/beam, respectively, for the 3-, 2- and 
0.88-mm continuum maps; 0.56, 0.75 and 1.21 Jykms~^/beam, respectively, for the 4 — )■ 3, 
5—^4, and 7—7-6 CO velocity-integrated maps. Offsets are given with respect to the 
source position as determined by Herschel of RA 14:24:13.9, Dec. +02:23:04 (J2000). The 
white cross indicates the position of the source as determined from the 0.88-mm map and 
corresponds to RA 14:24:13.98, Dec. +02:23:03.45 (J2000). 
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Fig. 4. — Spectral energy distribution of ID 141 (all measurements from this paper: see 
Tabled]). The squares indicate the ground-based observations and the triangles the Herschel 
SPIRE data. For comparison, the infrared SED of the starburst M 82 is displayed, redshifted 
to 2; = 4.24 and normalized to the flux density of ID 141 at the observed wavelength of 500 /xm: 
the crosses show the photometric data taken from th e literature and the continuous line is 
the ISO LWS spectrum of M 82 (IColbert et al. 19991 ). The dashed line displays the best fit 
to the photometric data of ID 141 using a single component, optically thick dust model (as 
described in text). 
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Fig. 5. — Left Panel - Observed CO fluxes vs. rotational quantum number (CO line SED) 
for ID 141. The single- component LVG model (described in the text) corresponds to n{B.2) = 
^q3.9 p^-3 ^^^ y^.^ ^ 4Q j^ j^^g^^ p^^^i _ Com parison of the CO line SEDs of ID 141 and 



selected high-redsh ift SMGs: GN20 {z = 4.0 51 (ICariUi et al. 20101) . SMMJ16359 {z = 2.5) 
JWeifi et al. 2005b[ ) and AzTEC-3 {z = 5.3) jRiechers et al. 201oh . The CO line SEDs are 
normalized by their CO (1-0) flux density. 
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Fig. 6. — The ratio Lctt/-^ftr a-s a function of the far-infrared 
galaxies (JMalhotra et aL 200 ll: IColbert et al 19991 lUneer et al. 20001 
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Spinoglio et al. 2005 



Carral et al. 19941 1 and ULIRGs (jLuhman et al. 19981: iLuhman et al. 20031) in the local uni- 
verse, and starburst galaxies and AGN at 1 < 2; < 2 (IStacey et al. 2010l ) and z > 3 - see text 
for references. The source ID 141 is labelled (and underlined) as well as other well-known 
local starbursts (including M 82) and high-redshift sources. The source ID. 81 {z = 3.04) is a 
recent measurement based on SPIRE FTS observations (Valtchanov et al. 2011). For ID 141, 
the bar shows the range in far-infrared luminosity after correcting for an amplification factor 
10 < /iL < 30. 



